Geometric quantum gates with superconducting qubits 
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We suggest a scheme to implement a universal set of non-Abelian geometric transformations 
for a single logical qubit composed of three superconducting transmon qubits coupled to a single 
cavity. The scheme utilizes an adiabatic evolution in a rotating frame induced by the effective 
tripod Hamiltonian which is achieved by longitudinal driving of the transmons. The proposal is 
experimentally feasible with the current state of the art and could serve as a first proof of principle 
for geometric quantum computing. 
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I. INTRODUCTION 

Superconducting qubits, also known as artificial atoms, 
have emerged as a promising candidate to achieve quan- 
tum computing.^ The properties of these nanoscale sys- 
tems can be designed to a large extent, and systems have 
been found were many logical ope rati ons can be per- 
formed within the decoherence timeP^l Superconduct ing 
qubits can be coupled via thin-film microwave cavities 4 ^ 
to allow for two-qubit gates and ultimately for univer- 
sal quantum computing. Furthermore, they have the in- 
trinsic scalability of condensed matter systems and the 
high-precision measurement features of quantum optical 
systems. 

Error correction theory predicts that fault tolerant 
quantum computing requires of the order of 10 4 qua ntum 
operations with only a single error on average!^ Con- 
trary to classical computation, where gates and errors are 
discrete, in quantum computation many small errors can 
accumulate to an eventual bit or phase flip. Therefore, 
enormous accuracy for single gates is required. Typically, 
the control parameters cannot be controlled to such pre- 
cision and especially the exact timing of control pulses 
remains challenging. As a possible solution, holonomic 
quantum computing was proposedPln this case, the uni- 
tary transformations depend only on the path which the 
control parameters trace in parameter space, but not on 
their timing. Furthermore, random rapidly fluctuating 
deviations of the actu al path to the desired one cancel 
to the first orderPiHl Thus, the precision of holonomic 
quantum gates can possibly be considerably higher that 
the precision of the control parameters itself. 

Abelian holonomies, referred to as geometric phases or 
Berry phases, have been observed in a wide v ariety of 
systems including superconducting qubitsP^E^ The situ- 
ation is quite different for non-Abelian holonomies nec- 
essary for universal geometric quantum computing. De- 
spite several theoretical proposals^ 3 -*^, no such adiabatic 
transformation has been experimentally observed in su- 
perconducting qubits, nor in any other systems. Here, we 
present a scheme for the implementation of a non-Abelian 
holonomy which is feasible with the devices and methods 



used in current experiments on transmon qubitsP 

Our method is based on the much studied tripod 
Hamiltonian 



3 



H = aV(f2i|0)(i| +h.c.)=ti 



Q\ ^3 

000 
000 
000 



a) 



where the f2j(t) are the control parameters (usually re- 
ferred to as Rabi frequencies) and the matrix representa- 
tion is given in the basis {\i)}, i = 0, 1, 2, 3. The first pro- 
posal to observe non-Abelian transformations in trapped 
ions was based on Hamiltonian which is sufficient to 
implement an arbitray £/(2)-transformation s betw een the 
states |1) and \2) used as the logical qubit .1121321 Similar 
structures have been recovered in many quantum systems 
and similar implementations have been proposecpl" 2 ^, 
but yet without experimental verification. 

Based on recent experiments 2 ^ we propose a way to 
implement a universal set of single qubit non-Abelian 
geometric transformations in a system of three super- 
conducting transmon qubits coupled to the same cavity. 
Each transmon is composed of two superconducting is- 
lands connected by two Jos ephs on junctions, thus form- 
ing a superconducting loop.^i The control parameters 
are the magnetic fluxes through the loops of each trans- 
mon, which can be controlled individually allowing us 
to adiabatically drive the system along a control cycle. 
With realistic approximations we are able to obtain an 
effective tripod Hamiltonian in a rotating frame. 

The paper is organized as follows. In Sec. [TTJ we intro- 
duce the physical setup of the proposed experiment and 
in Sec. Ill we derive the effective tripod Hamiltonian. 
Section [IV] is devoted to numerical studies. We conclude 
the work in Sec. IVl 



II. PHYSICAL SYSTEM 

The superconducting qubits c onsid ered here are com- 
monly referred to as transmons Their structure is 
similar to a charge qubit, but they have a much larger 
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total capacitance Ce such that the ratio of the charging 
energy Eq — e 2 /(2Cz) over the Josephson energy Ej is 
much lower than unity. This results in a small charge 
dispersion of the energy eigenstates, which in turn leads 
to a significantly reduced sensitivity to charge noise and 
much longer decoherence times, typically of the order of a 
few microseconds. On the downside, they have a smaller 
anharmonicity compared with charge qubits. 

Our scheme includes three transmons with frequencies 
ei/2n, i — 1,2,3. The transmons are coupled to a cav- 
ity mode with frequency lo/2tt. The combined sys tem is 
described by the Tavis-Cummings HamiltoniarPEl 



H = HcrfdJa - 
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where gi/2n is the transmon-cavity coupling frequency, 
are the usual Pauli operators for the i-th transmon, 
and a and are the bosonic annihilation and creation 
operators for the cavity mode. To arrive at the above 
Hamiltonian, we used the rotating wave approximation 
(RWA) , assuming that the coupling strengths gi are small 
compared with the excitation energies, which will be the 
case throughout the paper. 

Furthermore, we neglected higher levels of the trans- 
mons which is very well justified, because we do not 
drive the transmons transversally and therefore do not 
induce excitations to higher energy levels, and we ini- 
tialize the system in the one-excitation subspace.^ In 
this case, the only states involved in the dynamics are 
{ I l^ffs); |O e Sff): |0<?e<7), \0gge)}, where the first states in 
the tensor product represents the photon number in the 
cavity mode and the states |e) and \g) are the excited 
and ground states of each transmon. The Hamiltonian 
restricted to this subspace, written in matrix form, reads 
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where Aj = Si — oj is the detuning of the i-th transmon 
from the cavity (see Fig. [IJ . 

The qubit frequencies £j, the detuning Aj, and the 
system-cavity coupling strengths gi depend on the 
Josephson energy Ej\(j>i) and thus on the controllable 
flux <f)i though the i-th qubit. The Josephson energy can 



be written as Ej 
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where Ef max is 



the maximum Josephson energy and 4>i i s in u nits of the 
flux quantum h/(2e). Explicitly, we hav d 2 | 3 | 26 J 
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where ki is a constant depending on the system parame- 
ters which can be determined experimentally. 

By changing the flux <f>i (t) through the i-th transmon, 
we can control the coupling strength gi(t) as well as the 
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FIG. 1: (Color online) The level structure of the lowest-energy 
states of the Hamiltonian Eq. ([2| are shown. The photon 
excitation number of the cavity is labeled by \n), while \g) and 
\e) are the ground and excited states of the three transmons, 
respectively. The couplings between the one-excitation levels 
gi are small compared to the detunings Aj. 



detuning Aj(i). We separate the the dominant constant 
contributions [denoted with superscript (0)] which de- 
fines the properties of the non-driven system from the 
small time-dependent ones, which are used to drive the 
system. We employ the notation 



Flux: <f>i{t) 
Detuning: Aj(£) 
Coupling: g^t) 
Hamiltonian: H(t) 
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We assume that the flux modulation is small compared 
to the flux quantum, i.e. 5<pi(t) <C 1, and use a first-order 
expansion in 5cf>i to obtain the time-dependent quantities. 
Because the cavity frequency is independent of the flux, 
we have 5 A, = fej. Together with gi(4>) oc y/ £i(^>) we ob- 
tain a useful relation between the coupling and detuning 
variations 



8gi(5<f>i) 
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which is valid up to first order in 5<f>i. Since typically 
gi <C Si, the driving via the flux has a much smaller 
effect on the transmon-cavity coupling gi than on the de- 
tuning A,, and therefore results mainly in longitudnal 
driving. However, the variations in the detuning induce 
transitions in higher order perturbation theory, which we 
refer to as indirect coupling. Whether, the direct or the 
indirect coupling gives the leading contribution to the ef- 
fective tripod Hamiltonian depends on the ratio A^ / 
and will be studied below in detail. 
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III. EFFECTIVE TRIPOD HAMILTONIAN 

In this section, we show how to obtain an effective tri- 
pod Hamiltonian from the Hamiltonian Eq. (3j). In par- 
ticular, the normal tripod approach which solely utilizes 
the driving of the off-diagonals of the Hamiltonian will 
not work for our situation, because our control over gi is 
rather limited. Nevertheless, we will find that driving the 
diagonals results in an indirect coupling of the different 
eigenstates of which is of the desired tripod form. 

To this end, we assume that the time dependent fluxes 
54>i(t) oscillate with the frequencies C0i/2ir and we write 

S(j)i(t) = Fi(t)cos[ojit + ifi(t)], 



5Ai{t) = Li(t)cos[uit + <pi{t)], 
8gi(t) = Ti{t) cos[u)it + <pi(t)], 



(6) 



where the adiabatically changing amplitudes Li(t) and 
Ti(t) are related with the externally controllable magni- 
tude of the flux oscillations Fi(t) by Eq. Q. To realize 
a universal set of single qubit transformations, also the 
relative phases fi{t) of the oscillations need to change 
adiabatically in timeP3 

Anticipating that SH(t) drives transitions between the 
eigenstates of H^-°\ we diagonalize H^-°\ Up to the first 
order in g^p / Af\ the eigenstates of in the basis 
{11999}, \0egg), \0geg), \0gge)} are given by 
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The Hamiltonian H(t) = i? (0) + 6H(t) in this basis assumes the form 
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Here, it is clear that Sgi(t) and 5Aj(t) have to oscillate with frequency = E^ to induce an effective coupling. Moving 
into the rotating frame with respect to the diagonal dominant contribution H dlag = diag{0, Ei,E 2 ,E 3 } and using Eq. 
^ we obtain 



H D (t) = h 
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Here, we defined the effective Rabi frequencies 
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where Eq. ([5]) was used in the second line. In the RWA we 
can drop the oscillating entries of Eq. ( 9 ) and we arrive 
at the desired tripod Hamiltonian Eq. ( 1 ) . 
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For negative detunings, i.e. uj > s\°\ the direct cou- 
pling due to Ti and the indirect coupling due to Li add up 
increasing the strength of the effective coupling. Depend- 
ing on the ratio between the detuning and the energy gap 
we have two different regimes. If |A| 0) | < ef ] we are in 
the small detuning regime and the second contribution 
dominates. If |A-°^| >• e- ^ we are in the large detuning 
regime and the first contribution dominates. Theoreti- 
cally, both regimes yield the tripod form of the effective 
Hamiltonian. Since the different regimes have different 
requirements on the experimental setup, which are read- 
ily available for small detuning regime, we study this in 
more detail in below. 

We used two approximations in the derivation of the 
tripod Hamiltonian. Firstly, gf^ <C Aj ' was needed to 
derive Eq. Although there exist higher order terms 
which might seem to destroy the ideal tripod structure, 
these can all be removed within the RWA. Neverthe- 
less, there are higher order terms resulting in an effective 
coupling slightly lower than suggested by Eq. ( 10 ) and 



an optimal driving frequency marginally different from 
Eq. jij. Secondly, the RWA requires L l < lu, w a\°\ 
Both relations limit the effective coupling strength of the 
indirect coupling, while the direct coupling is limited by 
g^ -C sf 1 ^ which was used to write down Eq. S. To 
demonstrate a holonomy with current experimental lim- 
itations (decoherence times, transmon-cavity couplings, 

flux driving), one needs to reduce the detuning A| to 
the edge of the validity of the above approximations. 
This is studied in the following section. 



IV. RESULTS 

In this section, we verify the validity of above analyt- 
ical results with numerical studies. As an example, we 
choose to work with a particular implementation of non- 
Abelian operator proposed in Refs. 119)201 in which the 
Rabi frequencies fi^ are real and parameterized as 



fix = f) sin /3 cos a, 
f^2 = Osin/3sina, 
SI3 = f) cos p. 



(11) 



Accordingly, in the driving fields in Eq. ^ we take ipi = 
0. We assume f) to be constant, while the angles a, and 
j3 change adiabatically in time. 

Let us discuss some general properties of the tripod 
Hamiltonian in Eq. ([!]). It has two non-degenerate 
eigenstates, usually refer red to as the bright states, 
with energies Ml — ±?j^/f2f + 0| + f2| and, more im- 
portantly, a degenerate zero-energy subspace £(t) — 
span{|Di(t)} , |£>2(t))}. These so-called dark states 
are |-Di(i)) =cos/3(cosa|l) + sina|2)) — sin/3 |3) and 
l-CM*)) = cosa |2) — sin a 1 1). The system state \ip(0)) is 
initially prepared to be in this subspace, and if the control 
parameters are changed adiabatically, then the system 




FIG. 2: (Color online) The parameters f3 (black) and a (blue) 
as a function of time during the adiabatic control cycle. 



will stay in this subspace during the evolution. In partic- 
ular, for a cyclic Hamiltonian H(0) — H(T) the system 
will return to the initial subspace £ (0). However, within 
this subspace the state will undergo a non-tri vial U {2) 
transformation which is the holonomic operator! 28 * 29 ! 

The evolution in the parameter space begins and ends 
at the point (fii, f)2, ^3) = (0,0, fi) and by writing ex- 
plicitly the dark states using Eq. (11), we obtain that 



the initial zero-energy subspace is spanned by {|1), |2)}. 
These states are used as basis states of a logical qubit. 
An adiabatic change of the angles according to 

(a{t),0(t)) : (0,0) -> (0,7r/2) -► (7r/2,7r/2) -> (tt/2,0), 

(12) 

results up to a phase factor in a holonomic NOT gate 
[7 h oi = |1) (2| - 1 2) (1| for the logical qubit. Note that be- 
cause of the spherical parameterization, the Hamiltonian 
is cyclic. For better adiabaticity, we change the angles 
smoothly using sine functions and constants as shown in 
Fig.0 " 

To implement this loop in our setup, one can control 
the flux driving amplitudes Fi(t) and hence the longitu- 
dinal driving amplitude Lj (t) which is directly related to 
the Rabi frequencies by Eq. (10). The reference basis 
is given by Eq. Q and the initial degenerate subspace 
corresponds to |lj==«i and \2) = V2- In particular, we as- 
sume an initial state |1) and, for an ideal transformation, 
the final state is simply \ipidc&i{T)) = |2). Thus, the fi- 
delity, defined as F(t) = | (V'idcai (*) I V'C*)} | 2 > after the gate 
time T is the population of 1 2) . For clarity, we only show 
the populations of |1) and \2) in the figures below. 



We integrate the dynamics of the system using the 
ideal tripod Hamiltonian in Eq. ([I]) for different gate 
times T, as shown in Fig. (3^a) and (b). The fidelity is 
plotted over the gate time T in Fig. J^c) and shows the 
expected approach to unity in the adiabatic limit. The 
slightly oscillatory behavior observed in Fig. [3^c) is typ- 
ical for adiabatic gatesP^ However, one should not rely 
on local maxima of this curve, as their position depends 
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FIG. 3: (Color online) The population of the states 1 1} (black) 
and 1 2} (blue), obtained from the numerical integration of 
the Schrodinger equation using the effective Hamiltonian in 
Eq. (if with Q./2-k = 10.5 MHz. (a) The gate time T = 0.5 ^s 
shows good population transfer, while (b) for T = 0.3 g,s 
the smaller fidelity indicates the invalidity of the adiabatic 
approximation. Panel (c) shows the gate fidelity as a function 
of the gate time T. 



on the precise value of several experimental parameters. 
Instead, one should use gate times long enough such that 
even a local minimum provides a sufficiently good fidelity. 

In Fig. |4j we integrate the exact Hamiltonian of 
Eq. ([3]). By comparing the results with Fig. |3j we can 
judge whether approximations such as the rotating wave 
approximation are satisfied. We use parameters which 
yield an effective coupling of Q/2ir = 10.5 MHz (10 MHz 
indirect coupling and 0.5 MHz direct coupling). The re- 
sults follow closely to the ones obtained by the effec- 
tive tripod Hamiltonian, i.e., a gate time of 0.5 fis re- 
sults in a fidelity of almost unity, whereas a gate time of 
0.3 fis is not enough to justify the adiabatic approxima- 
tion. The fidelity is plotted as a function of the gate time 
in Fig.Qc), which shows much resemblance to the corre- 
sponding Fig. [3^c) except for a slight rescaling of the gate 
time. The reason for this rescaling is that /A^ = 0.2 
is not small enough to perfectly justify the ap prox imation 



g^ <C A^ u; , and therefore the formula Eq. (10 1 slightly 



,(0) 



overestimates the effective coupling as described in the 
previous section. 

Because the decoherence time of transmons to date is 
of the order of a microsecond, one would wish to increase 
the effective coupling to achieve faster holonomies. This 
can be done in various ways as shown in Fig. [5ja)-(d), 
always such that Eq. ( 10 ) suggests roughly double the 
effective coupling compared with Fig. |4]jc). One would 
expect good fidelities in half the gate time compared with 
Fig. [4](c). The easiest and readily available way is to de- 
crease the detunings A$ as shown in Fig. [5ja). However, 
the gate fidelity is by no means as good as expected from 



FIG. 4: (Color online) Panels (a) and (b): The population 
of the states |1) (black) and |2) (blue), obtained from the 
numerical integration of the Schrodinger equation using the 
exact Hamiltonian Eq. d3j. The parameters are as follows: 



u/2tt = 5 GHz, g{ 0) /2n = 60 MHz, g$" /2tt 



(°). 



MHz, 



gl '/2ir = 100 MHz, A K "> /2x = -300 MHz, A^> /2-n = -400 
MHz, A ( 3 0) /2ty = -500 MHz. The longitudinal driving has a 
restriction of max(Li) = 100 MHz, resulting in an effective 
coupling Q/2n = 10.5 MHz. The gate times are T = 0.5 /xs 
in panel (a) and T = 0.3 fj,s in panel (b). Panel (c) shows the 
fidelity as a function of the gate time T. 



the effective coupling. The reason is that the conditions 
Li <C A,- ' and gi <C A^ ' 1 used in the derivation of the 
effective Hamiltonian tend to get violated for decreasing 
A' - 1 . The situation is slightly better in Fig. pi b) and (c), 
where the higher effective coupling is achieved by increas- 
ing the transmon-cavity coupling gf^ and the driving 
amplitude L,, respectively. The only way to increase the 
effective coupling without affecting the validity of the ap- 
proximations is to simultaneously increase g\ , L^, and 
A[ \ which is shown in Fig. [^d). However, it may be 

hard to achieve such high transmon-cavity couplings g^ 
and driving amplitudes Li with current setups. 

We would like to add a note on experimental feasi- 
bility. The parameters used in F ig. [4] are realistic in 
an existing experimental setup^^ and can be achieved 
[see Eq. Q], for example, by using transmons with 
charging energy Eq — 2ttH x 280 MHz, Josephson en- 

6 (0) 



ergy E Jn 



2nh x 224 GHz, fluxes 



0.48426, 



JO) 



0.48489, and 



4 0) 



0.48550, and variations of 
the fluxes of up to max(Scj)i) — 6.3 x 10~ 4 . To verify 
the geometric transformation, one has to be able to read 
out the final state of the system. For this purpose, one 
increases the detunings considerably such that the en- 
ergy eigenstates Vi are approximate product states of the 
cavity and the three transmons. Then, it is sufficient to 
perform state tomography of the first and second trans- 
mon because the holonomy is a transformation between 
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qubits, which could serve as a first step towards geomet- 
ric quantum computing. Although we did not explicitely 
take into account the decoherence in our studies, the gate 
time is within the decoherence time for current experi- 
mental setups allowing for proof of principle experiments. 
The detailed effects of decoherence can be studied along 
the lines of Refs. l3~4"H3"o1 and will be presented in a future 
publication. We note that there could be considerable 
technical improvements in the near future concerning the 
decoherence time as well as the driving strength, leading 
to the possibility to carry out extensive small-scale quan- 
tum computing. We used the NOT gate as an example 
to calculate the gate fidelity, but the proposed scheme 
can_be utilized to carry out any single-qubit transforma- 
tio: 



FIG. 5: (Color online) Gate fidelity as a function of gate 
time T. Parameters are choses as in Fig. [4] but with stronger 
effective coupling achieved by using half the detuning Af (a), 
double the transmon-cavity coupling gf ' (b), double the driv- 
ing Li (c), and double of all of them (d). This results in the ef- 
fective couplings Q./2-k = 20.5 (a), 21 (6), 21 (c), 22 (d) MHz. 



v\ and i>2 only. State tomography has been demonstrated 
for up to three transmons in Refs. I32|33l 



V. CONCLUSION 

We proposed an experimental scheme for geomet- 
ric non-Abelian single-qubit gates with superconducting 
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